Nonlocality effects on Color Spin Locking condensates 
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We consider the color spin locking (CSL) phase of two-flavor quark matter at zero temperature for 
nonlocal instantaneous, separable interactions. We employ a Lorentzian-type form factor allowing 
a parametric interpolation between the sharp (Nambu-Jona-Lasinio (NJL) model) and very smooth 
(e.g. Gaussian) cut-off models for systematic studies of the influence on the CSL condensate the 
deviation from the NJL model entails. This smoothing of the NJL model form factor shows ad- 
vantageous features for the phenomenology of compact stars: (i) a lowering of the critical chemical 
potential for the onset of the chiral phase transition as a prerequisite for stability of hybrid stars 
with extended quark matter cores and (ii) a reduction of the smallest pairing gap to the order of 
100 keV, being in the range of values interesting for phenomenological studies of hybrid star cooling 
evolution. 
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I. INTRODUCTION 



Recently, the investigation of color superconducting 
phases in cold dense quark matter has received much at- 
tention [l|, 0, y, 0, [a, [a] J in particular due to thepossible 
consequences for the physics of compact stars [3,|g]. From 
the point of view of observational constraints on quark 
matter and color superconductivity in compact stars the 
cooling characteristics play a central role. It has been 
shown that the occurrence of a normal quark matter core 
would lead to a conflict with observations since the di- 
rect Urea (DU) process in normal quark matter would 
lead to enhanced cooling in disagreement with the data 
[3, [i3|- The DU conflict would be solved provided no 
ungapped quark modes occur in the quark core. This 
has been demonstrated on the example of a hypothet- 
ical pairing channel (X-gap) for the quark color which 
is ungapped in the 2SC phase (2SC+X phase) [13, [HI- 
However, the microscopic origin of the X-gap could not 
yet be specified. A microscopically well-defined pairing 
pattern which could solve the quark DU problem would 
be the CSL phase |12| corresponding to a spin-one con- 
densate [l3l . uA. [l5j . A prerequisite for the realization of 
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this pairing pattern in quark matter would be a sufficient 
flavor asymmetry to prevent the u-d pairing in the other- 
wise dominant scalar diquark channel of the 2SC phase. 
It has been demonstrated that under neutron star condi- 
tions the 2SC phase is indeed rather fragile and may not 
be realized for moderate coupling strengths [la|. Thus 
the CSL phase becomes particularly interesting for the 
solution of the quark DU cooling problem, and corre- 
sponding simulations will be performed as soon as the 
cooling regulators such as emissivities, specific heat and 
thermal conductivity will be provided. First steps in this 
direction have been made recently [o, ^M, ■ 

Most of the calculations of QCD superconducting 
phases have been done using the sharp cut-off NJL model 
(see Ref.[l[ and references therein). However, lattice 
QCD calculations [l0| indicate that quark interactions 
should act over a certain range in the momentum space, 
and various approaches to include nonlocality effects be- 
yond the NJL model have been suggested [20| . We refer 
to nonlocal separable interaction models as introduced, 
e.g., in the works [U, [H, [H, [H, [H, ^ and references 
therein, where it has been concluded that smoothing the 
cutoff leads to a reduction of the chiral condensate and a 
lowering of the critical temperature for the chiral phase 
transition. The question arises for the effects of nonlocal- 
ity on the spin-one gaps, to be explored by varying the 
form factor of the quark interaction from a sharp cutoff 
in the NJL model to smoothly decreasing form such as a 
Gaussian. First exploratory calculations reported in [23| 
have shown that the nonlocality could lead to a sizeable 



decrease of the energy gaps. 

In this paper, we investigate the robustness of CSL 
pairing against a modification of the sharp cut-off (NJL) 
in a systematic way by employing a separable, instan- 
taneous interaction with a Lorentzian-type interaction 
which allows to interpolate between the NJL case and 
very smooth interaction form factors of, e.g., the Gaus- 
sian type. This investigation is performed on the basis 
of recently developed parameterizations for the instanta- 
neous three flavor case u3 . 



II. NONLOCAL CHIRAL QUARK MODEL FOR 
THE COLOR-SPIN LOCKING (CSL) PHASE 

We investigate a nonlocal chiral quark model in which 
the quark interaction is represented in a separable way by 
introducing form factor functions g(p) in the bilincars of 
the current-current interaction terms in the Lagrangian 
[la . l2ll . l29j . It is assumed that this four-fermion inter- 
action is instantaneous and therefore the form factors 
do not depend on the energy but only on the modulus 
of the three momentum p ~ \p\. The ansatz for the s- 
wavc, single flavor diquark condensate characterizing the 
CSL phase as introduced in Rcf. [I4I ^ is a scalar prod- 
uct (locking) of the three-vector of antisymmetric color 
matrices (A2, A5, Ay) with the three- vector of Dirac spin 
matrices (73, 72, 71). Thus, the corresponding gap matrix 
A for the CSL phase reads 



A(73A2 -I-72A5 -f 71A7) 



(1) 



Since the two flavor channels decouple, the quark thermo- 
dynamical potential can be decomposed into single-flavor 
components 



n,{T,{^,f}) = J2 ^iT,fij) 



(2) 



f=u,d 



and it is sufficient to consider in the following the con- 
tribution of a single flavor only, which in the mean field 
approximation is given by 



n{T,^i) 



SG SH^ 



^E 



d^p 1 

(27r)3 2 



Trln 



i^-S ^{iuJn,p) 



.(3) 



where ^ stands for the chemical potential of that flavor. 
The first two terms are quadratic contributions of the 
mean field values 4> and A of the order parameter fields 
that signal chiral symmetry breaking and CSL supercon- 
ductivity, respectively. Their denominators contain the 



^ Note that such ansatz differs from the one in [l5l |. where the 
states are constructed with total angular momentum one. 



coupling constants G and Hy in the corresponding chan- 
nel. In ([3]) the sum is over fermionic Matsubara frequen- 
cies uin ~ (2n + 1)ttT, and the trace is over Dirac, color 
and Nambu-Gorkov indices. 

In our nonlocal extension, the inverse fermion prop- 
agator differs from the NJL model case by momentum 
dependent form factors g{p) modifying the mesonic and 
diquark mean fields 



s-\p) 



M7° - M{p) g{p)A 



where M(p) is the dynamical quark mass function 
M{p) = m + .g(p)(/) . 



(4) 



(5) 



Note that although the first two terms in Eq. fS]) do 
not have any explicit dependence on the form factors, the 
quantities </> and A do depend implicitly on them through 
the gap equations, see Eq. p^ below. 

After evaluation of the trace [l^l and Matsubara sum- 
mation the thermodynamical potential takes the form 



n{T,ti) 



(6) 

where E^ [p) denote the excitation energies for the modes 
k = 1 ... 6. The odd (even) indices denote particle (an- 
tiparticle) excitations corresponding each to a triplet of 
spin-one eigenstates. All modes have a gap in the excita- 
tion spectrum and can be brought into a standard form, 
which for Ei [p) reads 



8G 8Hy 






Ekip) + 2T\n{l + e-'^-^P'^/'^) 



El{p) = (eeff (P) - A*cff(p))' + A2ff(p) 
with the effective quantities 



ecff(p) = ^Jp'^ + M^sip) : 



Meff(p) 



Moff(p) 



M{p) , 



Mcff(p) = AiVl + AVW^" , 
Aeff(p) = ^^Ag{p) , 

and for E^^^iji) is given by 

El^ip) - {e{p)-pf + a^,^{p)A^g^{p) 
with the momentum-dependent coefficients 



(7) 

(8) 
(9) 

(10) 

(11) 

(12) 



03,5 (p) 




p^ Y 



e2(p) 



(13) 



where e{p) = y/p^ + IvP (p) . The remaining modes 
E2A,6ip) ^^^ obtained from £'1,3,5 (p) by changing ^ — > 
—fi in Eqs. ((7)l- p^ . Note that the modes £1,2 (p) corre- 
spond to the vanishing z-projection of the spin, Sz = 0, 
thus being inert against an external B-field. The remain- 
ing modes corresponding to 5*^ = ±1 are expected to get 
shifted (Zeeman effect). 

For given values of T and /i, the global minimum of 
ri(T, yu) in the space of the order parameters (j) and A 
corresponds to the thermodynamical equilibrium state. 
We obtain this state by comparing solutions of the gap 
equations 



S6 ~ ^A 



= 



(14) 



We present results for the case of vanishing temperature 
and finite chemical potential in the next section. 



III. MODEL CALCULATIONS 
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A. Form factors and their parameters 

In (j4]) and ([5]) we have introduced the same form fac- 
tors g{p) to represent the nonlocality of the interaction 
in the meson {qq) and diquark (qq) channels. In our cal- 
culations we use the sharp cutoff (NJL). Lorentzian with 
integer parameter a (La) and Gaussian (G), form factors 
defined as 



ffNJL(p) 

9ha{p) 
9g{p) 



6(1 -p/ A), 
[l + {p/Af-Y 
exp(-p^/A^) 



a > 2 



(15) 
(16) 
(17) 



where A is a cut-off parameter. These form factors are 
plotted in Fig. [TJ We achieve deviations from the NJL 
case (step function) by using the Lorentzian form factor 
with decreasing the a parameter. The Gaussian form 
factor appears on the other limit having a very soft mo- 
mentum dependence. 

To perform numerical calculations one has to specify, 
for each form factor, the following set of parameters: the 
light quark current mass (m), the coupling strength (G) 
and the range of the interaction (A). The diquark cou- 
pling constant H^ is fixed to the ratio H^/G = 8/3 in ac- 
cordance with the result of the Fierz transformation for a 
one-gluon exchange interaction. In this work we use the 
parameterizations recently given in Ref. [28[ and listed in 
Table |T1 They have been obtained by fitting the vacuum 
properties of the pion (/^ = 92.4 MeV, M^ = 135 MeV) 
and the vacuum constituent quark mass at zero momen- 
tum, AI = rn + 0. For the latter the phenomenologically 
reasonable values M = 330 and 400 MeV are been used. 

Note that the results to be presented below do not 
depend on the choice of the Lorentzian-type function as 
interpolating form factor. In fact, similar results have 
been obtained using other interpolating functions as, e.g., 
the Woods-Saxon form factors the parameterization of 
which is given in [28[ . 



FIG. 1: Form factors used to represent the nonlocality in the 
momentum space. Decreasing a in the Lorentzian-type form 
factor (notation La) causes deviations from the NJL model 
in a systematic way. 



M[MeV] 


Form factor 


A [MeV] 


GA^ 


m[MeV] 


330 


NJL 


629.5 


2.17 


5.28 




LIO 


649.2 


2.36 


4.71 




L5 


666.5 


2.49 


4.09 




L3 


685.8 


2.59 


3.25 




L2 


703.4 


2.58 


2.37 




G 


891.1 


3.88 


2.18 


400 


NJL 


587.9 


2.44 


5.58 




LIO 


600.3 


2.64 


5.01 




L5 


609.3 


2.78 


4.39 




L3 


616.2 


2.87 


3.55 




L2 


617.8 


2.83 


2.65 




G 


756.1 


4.22 


2.60 



TABLE I; Parameter sets for the nonlocal chiral quark model 
with Lorentzian and Gaussian form factors and for the NJL 
model. Sets for different fixed M = 330, 400 MeV are listed. 



B. Quark mass and CSL pairing gap 

First; we analyze form factors which do not deviate 
strongly from the NJL case, i.e. La for a > 3, shown 
as the grey area in Fig. [T] In Fig. [5] we compare the 
solutions obtained for the mass and the CSL gaps for 
two different sets of regularizations for fixed constituent 
mass: M = 330 MeV (left) and M = 400 MeV (right). 
For the parameterizations with a larger constituent mass 



in vacuum, one obtains a larger critical chemical potential 
/Xc for the phase transition from the chirally broken phase 
to the restored one, where the CSL pairing can occur. 
The gaps at the onset, A(^c), are larger whereas the 
mass gaps after the chiral transition are smaller. 



is also a characteristics of confi ning quark models within 
the Dyson-Schwinger approach [34| . Smoother form fac- 
tors could thus be more appropriate to model QCD inter- 
actions. Note, however, that the instantaneous nonlocal 
models with such smooth form factors lead to unrealis- 
tically large values of the chiral condensate (above 280 
MeV) in the vacuum [H]. 
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FIG. 2: The dependence of the dynamical mass Af (0) and the 
CSL pairing gap A on the chemical potential /j, for different 
NJL-like form factor models (NJL and Lorentzian with a > 
3). The parameterizations correspond to a fixed constituent 
mass: M = 330 MeV on the left and M = 400 MeV on the 
right panel. 



Like in the NJL case [12], the CSL gaps are strongly 
increasing functions of /i in the range that is relevant 

for compact stars, ^c < M ^ 500 MeV, where the upper 
limit is due to the threshold for the occurrence of strange 
quarks which allow pairing patterns like the CFL phase, 
more favorable than CSL (for recent phase diagrams for 
neutral matter in the three flavor case see [30, [3l|j [23 ) • 

On the other hand, Fig.[2|clearly shows that the pairing 
gaps in this nonlocal extension are reduced relative to the 
NJL ones: the smoother the form factor, the smaller the 
gap. The reduction could be up to a factor three in the 
case of Lorentzian with a = 3. 

But perhaps one of the most important effects of non- 
locality from the point of view of the phenomenology 
of compact stars is the shift of the chiral phase transi- 
tion to lower values of fi. A lowering of the critical den- 
sity for deconfinement makes stable hybrid star configu- 
rations with large quark matter cores possible [23, [331 ■ 
On the other hand, results from lattice QCD simulations 
for the quark propagator [19| show a very smooth four- 
momentum dependence of the quark self energies which 



FIG. 3: Same as Fig. 2 for smooth form factor models 
(Lorentzian with a = 2 and Gaussian). 



Therefore we have subdivided our discussion of dif- 
ferent form factors into two groups: those which lead to 
deviations from the NJL results within one order of mag- 
nitude and those resulting in larger deviations from the 
NJL model case (Lorentzian with a = 2 and Gaussian), 
see Fig. [3| For example, the shift in the critical chemi- 
cal potential for the onset of the chiral phase transition 
relative to the NJL case is less than 20 MeV within the 
first group, but larger 30 MeV for the second group, see 
Figs. [2| and m However, the qualitative behaviour of the 
chiral and CSL gaps is not affected by the choice of the 
form factors. 

To obtain the above results, we have kept fixed the 
ratio Hy/G at the standard value obtained from Fierz 
transforming the one-gluon exchange interaction. How- 
ever, in order to estimate the effect of possible uncertain- 
ties in this value, we have considered also the situation 
in which this ratio is taken to be twice its Fierz value. 
The corresponding results for 2Hy and all Lorentzian- 
type form factors under consideration are shown in Fig.|4| 
It is worth noticing that while the increase of Hy by a fac- 
tor 2 increases the CSL gaps by one (M = 400 MeV) or 
two (M = 330 MeV) orders of magnitude, the qualitative 



behaviour of M and A as functions of /x remains practi- 
cally unchanged. The lowering of /Xc and A as the form 
factor becomes smoother is also similar to the previous 
case. 
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FIG. 5: Excitation energies in the CSL phase as a function of 
the momentum for both parameterizations of the LIO model 
from Table I at the corresponding critical values of the chem- 
ical potential, fi = fic- Left panel: M=330 MeV, fic = 338 
MeV and right panel: M = 400 MeV, fic = 375 MeV. An- 
tiparticle modes E2, E4, and Ee are not shown. 



FIG. 4: Same as Fig. 2 and Fig. 3 but the coupling constant Hv 
in the diquark channel is doubled with respect to the usual 
value coming from Fierz transformed one gluon exchange in- 
teraction. 



C. Quasiparticle excitation spectrum 

In Fig.[5|wc show the quasiparticle excitation spectrum 
at the critical chemical potential for both parameteri- 
zations, respectively. We observe that the quasiparticle 
mode with the lowest energy band corresponds is Ei{p) 
with a minimum 



-El. mill = min[£;i(p)] 
p 



(18) 



being the most relevant quantity for possible applications 
of the CSL phase of quark matter to compact star cooling 
phenomenology. This minimum occurs at the Fermi mo- 
mentum p = pp. In a very good approximation pp can 
be represented by the lowest orders of a series expansion 
in the parameter s = ppg'pi^{pp)/gLa(j>F), which is a 
measure for the influence of the form factor 



pI 



MoffbF) 



-2A2g(p^) 



-Ois') 



M^siPp) 
M{pp)[M{pp) 



MUpf) 



NP[pp)-M^^{pp)l,i^ 



In the same order of the expansion in s, we obtain for 
the minimal excitation energy 



£;i,niiii = Aeff(pF)+0(s2) 



(20) 



(19) 



Although i?i,min rniglit be quite small (see below) it never 
vanishes. In fact, as in the NJL case |12l |. in the present 
class of models none of the dispersion relations lead to 
gapless modes. It is interesting to note that this low- 
est energy mode Ei (p) , relevant for compact star cooling 
phenomenology, is the one which is inert against the in- 
fluence of a strong external magnetic field typical for neu- 
tron stars since it belongs to vanishing spin projection, 
5. = 0. 

In Fig.|6]and Fig.|7]wc plot i?i,min as a function of /U for 
different models from the NJL-likc and the smooth form 
factor groups, respectively. The calculations are made for 
both sets of parameterizations of Table 1, corresponding 
to constituent masses of M = 330 MeV and M = 400 
MeV. According to the analytical approximative result 
of Eq. (|20|) . the behavior of the minimal excitation en- 
ergies can be understood as a product of the increasing 
/i-dependencc of the CSL pairing gaps and the decreas- 
ing one of the other factors in Eq. (fTT|) . For the param- 
eterizations with M ~ 400 MeV on the right panel we 
obtain that £^i,min is a decreasing function of n since the 
increase in A(/x) cannot overcompensate the decrease in 
M{pp)q{pp)/fj,cfr(pp). A similar effect has been reported 
for NJL models [i2|. On the other hand, for parameter- 
izations with M = 330 MeV the interplay between A(/x) 



and the other factors in Eq. pT|) is density dependent: a 
slightly increasing behavior of i?i.min at low densities is 
followed by a tendency to a saturation or even decreas- 
ing behaviour at high densities. For the group of NJL- 
likc form factors, our results for -Ei^min lie in the range 
of 50 — 500 keV and for the case of smooth form factors 
they are between 1 and 100 keV. 
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FIG. 6: Minimal excitation energies -Ei,min in the CSL phase 
as functions of the chemical potential /i for different NJL-like 
form factors. Parameterizations correspond to fixed M = 330 
MeV on the left and M = 400 MeV on the right. 



Our main results are summarized in Fig. [8l In the 
upper panel we show the scaling of the the critical chem- 
ical potential ^];°' for the onset of the CSL phase for the 
Lorentztian-type model La with the smoothness param- 
eter 1/a normalized to the NJL limit case, /i^'"^. In the 
lower panel we show a comparative plot of the minimal 
excitation energies i?]"nii,-, in units of the corresponding 
NJL counterpart Ef'^J^^^, evaluated at ^J^" and ^^''^, re- 
spectively. The corresponding results for 2Hy are also 
shown in Fig. [5] as open symbols. As we see the qualita- 
tive behaviour of both /i^^/^^jL ^^^ E\^^:^^/ Ef^J^:^^ as a 
function of 1/a remains unchanged. 

It is reinarkable that, as it would be expected from 
an expansion to lowest order in s, both quantities scale 
almost linearly with 1/a. In fact, a very good approxi- 
mation to our numerical results is obtained with 
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for a down to 2, where the slope parameters ^ and ^' do 
only moderately depend on the model parameterization 
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FIG. 7: Same as Fig.[6]for models that present a larger devia- 
tion from NJL models (Lorentzian with a = 2 and Gaussian). 



(Af) and the CSL coupling strength (Hy). For M = 330 
MeV we get ^ = 1.8 (1.5) and ^ = 0.13 (0.12), while 
for M = 400 MeV the corresponding values are ^ = 1.6 
(1.2), ^' = 0.17 (0.16). The numbers in parentheses are 
obtained by doubling Hy. 



IV. CONCLUSION 

We have studied the effect of instantaneous nonlocal 
interactions in the color spin locking (CSL) phase of 
quark matter. We have introduced momentuin depen- 
dent form factors to model the nonlocality and compared 
systematically with the local NJL counterpart. 

We have shown that there is a systematic lowering of 
the critical chemical potential for the onset of the CSL 
phase as well as for the minimal excitation energy (effec- 
tive CSL gap) as a function of the nonlocality which can 
be represented as a linear dependence on the smoothness 
parameter 1/a of the Lorentz-type form factor. These 
qualitative effects are shown to be robust under changes 
in the coupling constant used to represent the CSL inter- 
action. 

It has been found that hybrid star cooling requires all 
quark modes to be paired with a minimal pairing of the 
order of 10—100 keV to suppress the direct Urea pro- 
cess in quark matter. The present model for the CSL 
phase meets this requirement and calls for a more de- 
tailed analysis of the cooling phenomenology based on 
this microscopically justified pairing pattern. 

The smallest gap which governs the cooling phe- 
nomenology corresponds to the zero z-projection of the 
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FIG. 8: Two nonlocality effects on CSL condensates as a 
function of the smoothness parameter 1/a of Lorentzian-type 
models, for a — 2,3, 4, 5, 10, normalized to the NJL limit 
(a — > oo). First: the lowering of the critical chemical po- 
tentials fi]:°' respect to fj,^''^ (upper panel). Second: the 
reduction of the minimal excitation energies i?i',min respect 
to -Eilli^n (upper panel; values calculated at fi];°' and fi^''^, 
respectively). Results are shown as a function of 1/a for 
two different quark mass parameterizations in the vacuum: 
M = 330 MeV (triangles) and M = 400 MeV (squares). The 
open symbols show results for doubled CSL coupling strength 
(2 H,). 



spin and thus remains unaffected by the external mag- 



netic field of a compact star. Moreover, the CSL pairing 
pattern is a flavor singlet and insensitive to the flavor 
asymmetry in a compact star under /3-equilibrium. 

Therefore, the CSL phase with nonlocal instantaneous 
interactions is particularly interesting for applications in 
compact stars and allows to achieve a suitable description 
of quark matter properties by choosing the appropriate 
form factor models. Although it remains to be shown 
that hybrid star configurations with the CSL quark mat- 
ter phase could be stable, due to the small gaps, wc 
expect to have results reproducing those of the normal 
quark matter case, where stable quark matter cores in 
nonlocal models have been found. 
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